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ABSTRACT 


CODED THROUGHPUT PERFORMANCE SIMULATIONS 

FOR THE 

TIME- VARYING SATELLITE CHANNEL 
BY 

LI HAN, B.S. 


Master of Science in Electrical Engineering 
New Mexico State University 
Las Cruces, New Mexico, 1995 
Dr. William E. Ryan, Chair 


The design of a reliable satellite communication link involving the data transfer 
from a small, low-orbit satellite to a ground station, but through a geostationary 
satellite, was examined. In such a scenario, the received signal power to noise density 

I 

'/ 

ratio increases as the transmitting low-orbit satellite comes into view, and then 
decreases as it then departs, resulting in a short-duration, time-varying 
communication link. The optimal values of the small satellite antenna beamwidth, 
signaling rate, modulation scheme and the theoretical link throughput (in bits per day) 



have been determined in [1]. The goal of this thesis is to choose a practical coding 
scheme which maximizes the daily link throughput while satisfying a prescribed 
probability of error requirement. We examine the throughput of both fixed rate and 
variable rate concatenated forward error correction (FEC) coding schemes for the 
additive white Gaussian noise (AWGN) channel, and then examine the effect of radio 
frequency interference (RFI) on the best coding scheme among them. Interleaving is 
used to mitigate degradation due to RFI. It was found that the variable rate 
concatenated coding scheme could achieve 74 percent of the theoretical throughput, 
equivalent to 1.11 Gbits/day based on the cutoff rate R 0 . For comprasion, 87 percent 
is achievable for AWGN-only case. 
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Chapter 1 


INTRODUCTION 


We are concerned with reliable data transmission from a small, low-earth-orbit 
(LEO) satellite to a ground station, but through a geostationary satellite as in Fig. 1.1. 
The advantage of this approach is that a single ground station, which tracks only the 
geostationary (GEO) satellite, may be shared by a multiplicity of small satellites. In 
such a scenario, the received carrier power to noise density ratio, C/N 0 , will increase 
as the transmitting low-orbit satellite comes into view, and then decrease as it departs, 
resulting in a short-duration, time-varying satellite communication channel. C/N 0 and 
E/N 0 “profiles” are illustrated in Fig. 1.2, where E c is the energy per channel bit. The 
overall goal of the project is to design a modulation and coding scheme which 
maximizes the link throughput in bits per day for a given EJN 0 profile and a 
prescribed decoded error rate requirement. 

The whole project has been divided up into the following two sub-projects. The 
first sub-project [2] is a link design that covers the orbital simulation and the cutoff 
rate R 0 analysis, where the cutoff rate R 0 is a measure of the Shannon channel 
capacity from information theory. Because the simulation and analysis results from 
the first sub-project will be used in this thesis, we will discuss these results later in 
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this chapter. This thesis focuses on the second sub-project which is the error-control 
protocol design. The goal of the thesis is to find a near-optimal practical error-control 
coding scheme, where optimality is in the sense of maximizing the daily data 
throughput. 

In this thesis, the effect of high-level pulsed RFI is considered in addition to 
AWGN. In practice, a power-limited low-orbit small satellite might be transmitting to 
a ground station via a geostationary satellite repeater at the same time and using the 
same frequency band as intentional or unintentional transmissions from a ground- 
based emitter (e.g., radar). When the data bit duration is shorter than that of the RFI 
pulses, the degradation can be very severe. 

FEC coding schemes which concatenate convolutional with Reed-Solomon (RS) 
codes, with interleaving have already proved to be efficient counter-measures against 
RFI [3]. Both fixed and variable rate coding schemes will be examined in this thesis. 

1.1 Review of the First Sub-Project Results 

In the first sub-project, a LEO satellite with an orbital period of 102.86 minutes 
(14 orbits/day) and an inclination angle of 100 degrees was assumed. Further, the 
LEO satellite was assumed to have the following parameters. 

• Spin stabilization in a nadir orientation, i.e., the long axis of the satellite intersects 
the center of the earth. 

• Transmitter power: 3 W 
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• Antenna: circularly polarized, helical, non-gimbaled, pointing away from the 
earth’s center. And the following 3 dB beamwidths B are considered: 70, 52, 36, 
28, 20, and 14 degrees. 

• Frequency: S-band (2.2-2. 3 Ghz) 

The GEO satellite was taken to be NASA’s Tracking and Data Relay Satellite 
(TDRS) located at 41 degrees west, in geostationary orbit. The TDRS parameters 
assumed were: 

• G/T: 8.9 dB/K 

• Antenna: 16 foot parabolic, circularly polarized, and capable of open-loop 
tracking. 

The theoretical daily throughput, which is the number of information bits 
transmitted while E c (t)/N 0 > 0 dB* , based on the cutoff rate R 0 is defined as 

T R„ daily = JX(0 (bits) (1.1) 

one 

day 

where the integration is over the time in a day for which E c ( t)/N 0 > 0 dB. 

Assuming MPSK signaling on an AWGN channel with two-sided power spectral 
density Nq/2, the soft-decision expression for R 0 (t) is then given [4] by 


* In practice, carrier recovery for an MPSK signal is difficult below this value [5]. 
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' M-\ f r ( t \ f-n.lcW' 

. R*(0 = K ~ R Aog M £exp — rr^sin 2 — (bits/sec) (1.2) 

\t?o V N 0 v MJJ) 

E ( t ) 

where — — (E s = E c log 2 M) is the symbol SNR, and R c is the code bit rate. From 
N 0 

(1.1) and (1.2), we note that T Rn is a function of an E/Nq profile together with the 
signal set size M. 

_ The main result of the first sub-project shows that, for the given system scenario, 
the theoretical daily throughput T Ro Mly is maximized at 1.5 Gbits when the small 
satellite antenna beamwidth uses 20 degrees, the modulation scheme uses binary 
phase shift key (BPSK) or QPSK, and uses R c = 4.5 Mcbps as the code bit rate. We 
also note that, with 20 degrees antenna beamwidth, the LEO and GEO satellites have 
two contacts per day, which last about (3 min + 4.4 min) = 7.4 min while E c (t)/N 0 > 
OdB. 

In this thesis we estimate the throughput achievable by practical fixed and 
variable rate coding schemes, for both the AWGN and RFI situations. If we were to 
apply a code to the optimal ( B = 20 degrees, R c = 4.5 Mcbps, BPSK or QPSK) E/Nq 
profile of the first sub-project to simulate its throughput, the simulation time would be 
unmanageable because the number of code bits for this optimal case is approximately 
(7.4 min)(60 min/sec)(4.5 Mcbps) *2 x 10 10 code bits. Thus, we instead consider a 
single contact in the total C/N 0 profile whose form lends itself to analysis. 
Specifically, we shall assume a C/N 0 profile (actually, contact) of the form 
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(1.3) 


C(t)/N 0 = A q exp (-t 2 /2o 2 ) 
where the maximum value, A 0 , of C(t)/N 0 occurs at time t = 0 by our choice of time 
origin, and a is a measure of the duration of the communication link. Fig. 1.3 shows 
the Gaussian approximation of (1.3) for the C/N 0 profile. 

The C/No profile of (1.3) gives rise to an E/N 0 profile through the chosen code bit 
rate R c : 

E c (t)/N 0 = (Aq/R c ) exp(-t 2 /2o 2 ) (1.4) 

As mentioned, we shall assume communication is possible only when E c (t)/N 0 > 0 
dB. The interval [-r 0 , / 0 ] for which this is true is easily found by solving for / 0 in the 
equation E c (t)/N 0 = 1, yielding 


Actual C/No Profile vs. Gaussian Approximation 



Fig. 1.3 Gaussian Approximation of C(t)/N 0 Profile 
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t 0 = a ^/2 In (AJR C ) 


(1.5) 


which exists only if R c < A 0 . 

Then the theoretical throughput for this single contact is 

T Rb = )u*)dt (bits) (1.6) 

According to the analysis of [1], T Rb is maximized when R c = Aq/2 so that 
(W^o)Max = 4(Oy^ 0 «3dB. 

Hence, hereafter in this thesis we shall assume the optimal E c /N 0 profile is 

E c (t)/N 0 = 2 exp(-/ 2 /2a 2 ) (1.7) 

Further, we assume this profile is known a priori to the variable rate scheme since it 
is computable within reasonable accuracy from the known orbits of the 
communicating satellites. 

Ryan and Han [1] also showed that if t 0 corresponds to 15,000 8-bit RS symbols, 
the maximum T Ro at the optimal signaling rate R c and optimal modulation scheme, 
BPSK. or QPSK, was equal to 175,300 bits. This leads us to define the throughput 
efficiency S as g = T coded /T Ro = T coded /1 75, 300, where T coded is the coded 

throughput for a single contact. We shall use the value of t 0 indicated above and, 
hence, this expression for S through out the thesis. Note that once the efficiency for a 
particular coding scheme is estimated as above, it may be applied to the daily 
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theoretical throughput T Rn dally , to estimate the daily practical throughput. That is, the 
estimated daily throughput for code £ is given by 

daily ~ Zg T R„, daily 0 - 8 ) 

1.2 Outline of the Rest of Thesis 

The reminder of this thesis is outlined as follows. The system and channel model, 
including the RFI model, are described in Chapter 2. Then in Chapter 3, we 
numerically estimate the throughput efficiencies for the AWGN channel. Three 
coding schemes, one fixed rate and two variable rate schemes, will be examined. Both 
the numerical procedures and results will be discussed in detail. Chapter 4 will 
describe the effect of RFI on the best coding scheme of Chapter 3. We also discuss 
the convolutional interleaving technique. Finally, a summary and conclusions will be 
made in Chapter 5. In Appendix A, we will give a brief introduction to punctured 
convolutional codes, and in Appendix B we will describe all the programs used in this 
thesis, both in Matlab and in C, also give a sample simulation program, written in C, 
for the AWGN-plus-RFI channel. 
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Chapter 2 


SYSTEM AND RFI MODELS 


2.1 System Model 

The complete system transmission model is shown in Fig. 2.1. The left-hand side 
is the small, low-orbit satellite transmission system, the right-hand side is the ground 
station receiving system, and in between (e.g., AWGN plus RFI) is the physical time- 
varying channel. 

The FEC coding schemes employed in this system are concatenated codes which 
use a convolutional code as the inner code and an 8-bit-symbol RS code as the outer 
code. Because RFI will usually create a burst error, when this burst error covers more 
than one channel code bit, we add a convolutional interleaver between the inner 
convolutional code and the physical channel, which has the affect of scattering the 
RFI over widely separated code bits. The details of the convolutional interleaver will 
be discussed in Chapter 4. 

From the view of the outer RS code, the section between a-a and b-b can be 
considered as a super-channel. It is known that, the errors at the output of the Viterbi 
decoder occur in bursts [6], thus this super-channel is a burst error channel. We may 
improve the performance of the RS code by adding another (symbol) interleaver 
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Source Data 


Estimated Data 



AWGN RFI 


Fig. 2.1 System Model 
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between the inner code and the outer code to make the errors out of this super- 
channel appear as random as possible. 

The data transmission and reception process then proceeds as follows. The data at 
the transmitting terminal input are RS encoded, symbol interleaved, convolutionally 
encoded, and bit interleaved. The data are then BPSK (or QPSK) modulated and sent 
through the satellite channel. At the ground terminal receiver side, the received data 
need to go through the BPSK demodulator, bit de-interleaver, Viterbi decoder, 
symbol de-interleaver and RS decoder, in that sequence. We shall require that the 
probability that the RS decoder is unable to decode, P cw , be no more than 10‘ 5 : 

P < 10’ 5 

We shall call P cw the probability of code-word error for simplicity. 

2.2 Discrete-Time Equivalent Model for Simulation 

Because computer simulations are necessarily performed in discrete time, it is 
necessary to derive an equivalent discrete-time system model for our simulations. The 
next sub-section does this for the AWGN channel and the sub-section after that does 
it for the AWGN-plus-RFI channel. 

2.2.1 AWGN-Only Discrete-Time Model 

In order to derive the discrete-time model, we begin with a simplified BPSK 
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b k (0 or 1) 


A 



Fig. 2.2 Simplified System Model 

system model as in Fig. 2.2. Here b k is the k A information bit, CO 0 is the transmitting 
carrier frequency, T c is the code bit duration, n{t) is the AWGN continuous wave, i(t) 
is the pulsed continuous wave (CW) RFI, and c h n k and i k are the signal, noise and 

A 

RFI’ s values at sample times kT c Respectively. Finally, bk is the estimated data. 

When the additive noise is white and Gaussian with spectral density Nq/2, it is 
easy to derive [7] that n k is zero mean Gaussian with variance N 0 /2 and we write n k ~ 
r](0, N 0 /2). Further, it can easily be shown that c k = ± 1 so that E c = 1 . The resulting 
discrete-time AWGN-only channel model is then as shown in Fig. 2.3. 


A 


b k (0 or 1) 

Encoder 

c*-±l 

ft * 

Viterbi 

b k 


4 

J * 

k 

Decoder 









n k ~ tj(0, AV 2) 

Fig. 2.3 Discrete-Time AWGN-only Channel Model 
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2.2.2 AWGN plus RFI Discrete-Time Channel Model 


The RFI model can be fully characterized by RFI type, RFI arrival rate, RFI-to- 
code bit duration ratio P, and RFI power. In this section, we will first give the 
assumptions about these RFI parameters, then give the RFI model derivation, 
followed by some necessary RFI amplitude calibration. 

Assumptions 

• Only the pulsed CW RFI is considered in this thesis, and the CW RFI has the 
form 

i(t) = a ■ cos (© , t + (j) ) for x, <t<x 2 

= 0 otherwise (2.1) 

where a is the RFI’s pulse amplitude, co, is its radian frequency, <j) is its random 

phase, and x, is its starting time, x 2 is its ending time, thus, the RFI’s duration x is 

equal to x 2 - x,. 

• Two kinds of RFI arrival models, periodic and Poisson, are studied. For the 
periodic arrival RFI, a duty cycle is used to characterize its arrival rate. For 
example, if the duty cycle is equal to 10’ , the RFI starts after every 100 channel 
bits. The Poisson process, a very important counting stochastic process, is the 
other RFI arrival model. The definition of the Poisson process, denoted by N(t ) 
with parameter mean kt and arrival rate k, can be found in many books (such as 
[8], [9]). In this thesis, the thing we do need to know about Poisson process is not 
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the process itself but the time interval between two adjacent arrivals of RFI 
events, called the interarrival time. Fig. 2.4 shows the definition of the interarrival 
time, denoted by T n . Ochi [8] proved that the interarrival time of a Poisson 
process with arrival rate X obeys an exponential probability law with mean MX. 
Thus, the probability density function (pdf) of the interarrival time is 

p(T„)=\e-"- (2.2) 

Fig. 2.5 shows the example of both periodic and Poisson arrival models for RFI. 

• The RFI-to-code bit duration ratio (3 is defined as t IT C . Basically, this ratio 
tells us how many channel code bits will be affected by each RFI pulse. The 
total number of channel code bits affected by RFI is equal to either ["pi or Tpl + 
1 depending on the RFI’s starting time x,, where fpl is the ceiling faction of P 




time 

1 

event 


Fig. 2.4 Definition of Interarrival Time 


t, t 2 



Fig. 2.5 Examples of RFI arrival type (a) Periodic (b) Poisson 
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whose value leads to the next nearest integer from (3. For example (see Fig. 2.6), 
when P = 3.4, then RFI affects 

4 bits, if x, < 0.6 or 5 bits, if x, > 0.6 
In general, x, has the form X[ = K x -T c + K 2 , where K } is an integer and 0 < K 2 < 
T c . However in the above example, we have simplified the notation, letting X] = 

K 2 . 

• In the literature (e.g., [10], [1 1]), some authors assume that RFI has infinite power. 
We shall instead assume that the power of the RFI is 20 dB more than that of the 
channel code bit. 


Model Derivation 

Base on the above assumptions, (2.1), and focusing on the RFI part of the system 
model of Fig. 2.2 (see Fig. 2.7), the discrete-time RFI value i k can be derived as 

th • 

follows. When the interval (x ]5 x 2 ) does not include the k code bit, i k = 0. Otherwise, 


code bits: 


RFI (t, < 0.6): 


RFI (t, > 0.6): 


T r 



wL___ 



kU. 


! 

-U. ■ 

i 

i 




0.6 


ti 


L, 


P = 3.4, RFI affects 4 bits 


/t 


'i 

/ 


P = 3.4, RFI affects 5 bits 






Fig. 2.6 Examples of Total Number of Code Bits Affected by RFI 
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Fig. 2.7 RFI Model 


C iin(T 2 ,(A + l)7^.) | 2 

ci • cos(co /? + <())•/ ■— cos(co 0 1) ■ dt 

ax ( T t ,kT c ) 1/ / c 


(2.3) 


£ iin(T 2 ,( A: + l)7; ) a 
ax(t| y kT c ) 2 


• cos[(co , - (O 0 )•/ + (()]• dt 


r iin(T 2 ,(i + i)r ( .) a 
ax ( T , ,kT c ) 2 


(— -cos[(co 0 +(0 ,)-t+$]-dt 


The high frequency component (i.e., the second term) will be filtered out by the 
low pass filter (integrator), hence 


r iin(x 2 ,(*+])r c ) a 2 .. . 

r irT \ T'l Tr- C °s[(®/ -CD 0 ) -/ + <}>]- 

ax(T,,*r c ) 2 v T 


dt 


(2.4) 


Now defining Aw = (w ; - w 0 ) , and substituting into (2.4), we obtain 


2 ^ 'C03[4co 


(2.5) 


( 2 . 6 ) 


By assuming that the satellite system has a baseband bandwidth of 1/T C , we find it 
useful to write Aw as 
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Aco = 27i • ( a / T c ) 


for 0 < a <1 


(2.7) 


Finally, the discrete-time channel model for the AWGN-plus-RFI case is shown 
in Fig. 2.8. 

Calibration 

In this section, we calibrate the RFI amplitude, a, based on the system model we 
used and employing our assumption of an RFI-to-signal power ratio of 20 dB. We 
have 

P m - 100- P„„ K = m-(EJT c ) = 100/7; (since /:, = 1, (2.8) 

But P XF ,=~r-l'm 2 dl 

1 7 * 

= — | a 2 cos 2 (co , -t + §)dt 

T c 

^al 

2 

Substituting above into (2.8), we have 


A 


b k (0 or 1) 

Encoder 

c*=+l 
Uj 

p) 

A > 

Viterbi 

bk 


S 

/ s 

y 

Decoder 

> 


n k ~ r|(0, Ay 2) i k 

Fig. 2.8 Discrete-Time Channel Model 
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RFI 


a 2 _ 100 
2 “ T c 



(2.9) 


Now by substituting (2.7) and (2.9) into (2.6), we can further simplify the 
discrete-time RFI value i k to 



2 • 2n ■ a / T c 


•sin( 


2n a 

~T~ 


■t + 


♦>i 


min(t 2 , (A + l)7^) 
max(x , ,kT c ) 


_5_ 

an 



min(T 2 ,(& + 1)7 C ) +<j>]- sin[ 


2na 

T" 


max(x , , kT c ) + <j> ] 


( 2 . 10 ) 
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Chapter 3 


THROUGHPUT ESTIMATION METHODS 
AND RESULTS FOR 
THE AWGN-ONLY CHANNEL 


3.1 Selected Coding Schemes 

One fixed rate and two variable rate concatenated FEC coding schemes are 
examined for the AWGN channel in this thesis, as shown below: 

Fixed Rate Coding Scheme 

Code 1: rate 1/2, constraint length 7 convolutional (inner) code concatenated with 
a (255, 223) RS (outer) code over Galois Field (GF) (256) [6] (see Fig. 
3.1). This coding scheme is the NASA standard. 



RS 


Convolutional 



(255, 223) 


(k=7, r=l/2) 



Fig. 3.1 Concatenated FEC Coding Scheme 1 (Fixed Rate) 


Variable Rate Coding Schemes 

Code 2: variable rate (1/2, 2/3, 3/4, 4/5), constraint length 7 convolutional (inner) 
~ code concatenated with a (255, 223) RS (outer) code over GF(256) (see Fig. 
3.2). For convolutional coderates larger than 1/2, we assume punctured 
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, RS 


Convolutional (k=7) 

(255, 223) 


(i=l/2, 2/3, 3/4, 4/5) 


Fig. 3.2 Concatenated FEC Coding Scheme 2 (V ariable Rate) 


coding [12]. Also, we considered higher rate codes (e.g., 5/6, 6/7, ...), but 
the P cw was never achievable at these rates. 

Code 3 : variable rate (1/2, 2/3, 3/4, 4/5), constraint length 7 convolutional (inner) 
code concatenated with a (255, k= 223,225,. ..,255) RS (outer) code over 
GF(256) (see Fig. 3.3). 



Fig. 3.3 Concatenated FEC Coding Scheme 3 (Variable Rate) 


By using variable rate coding schemes, better throughput efficiency can be 
obtained, and hence better daily throughput. But the complexity of the system is 
increased. This tradeoff needs to be considered when a final decision which is made 
on which coding scheme is to be used. 

3.2 Throughput Estimation Procedures and Their Program Flow Charts 

Before we discuss the throughput estimation procedures, the relevant system 
conditions, requirements, and assumptions are summarized in Table 3.1. As 
mentioned, we will estimate the coded throughput efficiencies in one contact time 
interval [-/ 0 , t 0 ]. Since it is easier to perform our computations when time is in units of 


20 









symbol number, we let t 0 corresponds to 15,000 8-bit symbols. Then the E c /N 0 
profile becomes 

E c (s)/N 0 = 2 exp(-s 2 /2a) (3.1) 

where 5 is the symbol number, and 5 € [-15000, 15000]. The value of the variance can 
be determined from our earlier assumption of E c /N 0 = 0 dB at time s = ± 1 5000: 

2 exp(-15000 2 /2a 2 ) = 0 dB = 1 => a = 12,740 (3.2) 

Such an E c /N 0 profile is given in Fig. 3.4. 

Since the error rate performances of both the RS codes and the convolutional 
codes for the AWGN channel are well known (see below), it is possible to compute 
these quantities without actually simulating the system of Fig. 2.3. 


EcNo(dB) vs. symbol number for time-varying AWGN channel 



Fig. 3.4 Time-Varying Gaussian Shaped E c ( s)/N 0 Channel Profile 


* This number has to be large enough to accommodate many RS codewords. 
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Table 3.1 Summary of the Simulation System 


Items 

Requirements 

E c /N 0 Profile 

E c (s)/N 0 = 2 exp(-s 2 /2a 2 ), 
s e [-15000,15000], a =12,740 

Modulation Scheme 

BPSK 

FEC Coding Schemes 

Concatenated Codes 1 , 2, 3 

P cw Requirement 

P < 10' 5 

Assumptions 

E c /N 0 >0, Symbol Interleaver Used 


As indicated in Table 3.1, we will require a decoded outer RS codeword error 
probability of P cw < 10‘ 5 . Hence, the coded throughput T e is actually the total 
number of data bits corresponding to P cw < 10’ 5 . For code 3, the rate of the outer 
code and then the inner code was varied until P cw < 10' 5 was attained. 

The codeword error rate P cw is given [6] by 

255 

P„= ’ ( 3 - 3 ) 

j = e + 1 


where e = 


255 -k 
2 


is the error correction bound for the outer RS code. The quantity 


P is the 8-bit symbol error probability at the output of the convolutional decoder 
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which, assuming interleaving between the two codes, is related to the convolutional 
decoder bit error rate pasP = (\-p).ln turn, when Viterbi decoding is used, p is 
upper bounded [6] as 



where the convolutional code has rate k’/n’, minimum free distance dj, and 
information weight spectrum {b d }. Although E c is time-varying, we assume that it is 
relatively constant over the duration of a codeword as will be the case in practice 
where several thousands codewords are transmitted per second. 

We now consider the throughputs T a , a e {1,2,3}, of above three coding schemes. 
For both fixed and variable rate coding schemes, T a can be computed as follows 

7>8-(255)f>(c) a e {1,2,3} (3.5) 

C— ] 

where N cw is the number of codewords transmitted in the interval [-15000, 15000 ] 
which achieve P cw < 10" 5 and r(c) is the composite (convolutional plus RS) rate of the 

1 223 

c h codeword. For the fixed rate code, r(c) = for all c, so that T, = 

2 255 

(30,000)-(l/2)-(223/255) = 104,941 bits 

Finally, the throughput estimation procedure then consists of the following three 
steps: 

Step! : Generate a time- varying Gaussian shaped E c ( s)/N 0 channel as in Fig. 3.4, 
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where s is the symbol number in the interval [-15000, 15000]. 

Step 2: For each E c ( s)/N 0 value, calculate p by (3.4) and then calculate P cw by 
(3.3). If P cw < 1 O' 5 , increment N cw by 1. Otherwise, for code 2 and 3, adjust 
the rate(s) so that P cw < 10" 5 , and store r(c) for later use. 

Step 3; When all of the 30,000 symbols are exhausted* , calculate T a by (3.5), 
then compare them with the theoretical throughput T R = 175,300 bits to 
obtain their efficiencies. 

Fig. 3.5 and Fig. 3.6 show the flowcharts of the Matlab simulation programs for 
code 1 and code 3, respectively. 

In Fig. 3.5, we first set the maximum symbol value N^ mMox = 15,000, and 
initialize the first symbol N^JX) = -^symMax- £ c (s)/jV 0 is then obtained from (3.1). P cw 
is then calculated from (3.3) and (3.4) and compared with our threshold 10‘ 5 . Then we 
move to the next symbol and update the parameters for that symbol, and repeat the 
process until N^Jp) extends Finally, we calculate the coded efficiency as S a 

= T a /1 75,300, and save and plot the results. 

For code 3, the procedure is similar to the that of above, except we need to adjust 
(using higher rate) the code rates so that P cw is just less than 10" 5 . To do this, At the 
beginning, we put the convolutional code rates into a vector, denoted by code_rate lx4 
= [1/2, 2/3, 3/4, 4/5], and i is the index of code_rate lx4 . Then for each E c ( s)/N 0 


* The last codeword will typically extend beyond the last (30,000 th ) symbol 
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Fig. 3.5 Flowchart of Matlab Program for Fixed Rate Coding Scheme 
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Fig. 3.6 Flowchart of Matlab Program for Variable Rate Coding Scheme 
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value, we start using rate 4/5 first (i.e., start at i = 4), and set the code rate of RS K(c') 
= 223. Then we calculate p and P cw as before. If P cw < 10' 5 , keep the convolutional 
code rate but increment the RS rate, denoted by ( N-255 , K(c)), K(c ) by 2, and re- 
calculate P cw until K(c) > N or P cw > 10‘ 5 . We then store the current convolutional 
code rate and previous RS rate (i.e., K(c)-2) for 5th symbol (i.e. /'th codeword). If P cw 
> 10~ 5 , decrement convolutional code rate index i by 1 (i.e., move to the next lower 
rate), re-calculate P cw until P cw < 10" 5 or i=l, and then adjust K(c) as above and store 
the final code rates for the later throughput calculation. 

3.3 Numerical Results 

Fig. 3.7 plots the P cw profile for coding scheme 1, the fixed rate scheme. From 
the P cw profile, we note that P cw < 10* 5 for all 30,000 symbols, so that the throughput 
is 7j = (30,000)-(l/2) (223/255) = 104,941 bits, using the cutoff rate criterion, we 
observe the efficiency S\ =T\/ T Ru = 0.60. 

The P cw profile also shows us how powerful this NASA standard coding scheme 
is, so that we do not need to go further to check the performance of longer constraint 
length codes. But on the other hand, we note that P cw is much less than the prescribed 
threshold of 10°, which means that we waste a lot of coding overhead. In order to 
transfer as many information bits as possible, we try to use higher code rates in the 
higher SNR region, giving rise to coding schemes 2 and 3. 
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EcNo(dB) vs. symbol number for AWGN channel 




symbol number x 


Fig. 3.7 P cw Profile of Coding Scheme 1 for AWGN Channel 
Fig. 3.8 and Fig 3.9 plot the P cw profiles and the varying code rates for coding 
scheme 2 and scheme 3, respectively. As can be seen, the P cw of coding scheme 2 is 
getting close to the threshold , and the P cw of coding scheme 3 is almost right below 
the threshold line. The throughputs for these two codes computed from (3.5) yield T 2 
=144,740 bits, T 3 = 153,320 bits. The efficiencies are then S 2 — T 2 /T Rg = 0.83, S 3 = 

y T k = 0.87. This significant improvement is not surprising since the higher rate 
codes are used in the higher SNR region. 


28 





EcNo vs. symbol number for AWGN channel 




Fig. 3.8 P cw Profile of Coding Scheme 2 for AWGN Channel 
From these two figures, the varying of the code rate can be seen very clearly. In 
Fig. 3.8, we note that the convolutional code rate increasses when the E c /N 0 level 
increases. Further, in Fig 3.9, the RS rate may increase with E c /N 0 while the 
convolutional code rate remains constant. Hence, the composite code rate curve of 
coding scheme 3 is closer to the cutoff rate R 0 curve than that of coding scheme 2. 
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Table 3.2 summarizes the simulation results of the AWGN channel for the 


selected code schemes. As mentioned before, the throughput efficiency is obtained by 
considering a single contact of a total of 30,000 symbols, and the daily coded 
throughput is then obtained by multiplying the efficiency with the daily theoretical 
throughput of T R = 1.5 Gbits as determined in [2]. 


EcNo vs. symbol number for AWGN channel 



symbol number 

code rates and Ro bound vs. symbol number for AWGN channel 
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Fig. 3.9 P cw Profile of Coding Scheme 3 for AWGN Channel 
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Table 3.2 Summary of Simulation Results for AWGN Channel 


Coding 

Schemes 

Throughput 

Efficiencies 

Daily Coded Throughputs 
(Daily Theoretical Throughput = 1.5 Gbits) 

Code 1 

0.60 

0.9 Gbits 

Code 2 

0.83 

1 .2 Gbits 

' Code 3 

0.87 

1 .3 Gbits 
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Chapter 4 


EFFECT OF RFI ON THE BEST 
CODING SCHEME 


4.1 Simulation Procedure 

- Code 3, the best coding scheme from the previous chapter, will be examined first 
for the AWGN-plus-RFI channel to observe how severely the RFI can affect the 
throughput performance. Based on this result, it will be decided whether the 
throughput performance for Code 1 and 2 on AWGN-plus-RFI channel needs to be 
checked. 

Weinberg, [10] and [11], derived some formulas for the coded BER performance 
on the RFI channel, but unfortunately, these formulas are not easier than simulations. 
Thus we decided to obtain the coded throughput of AWGN-plus-RFI channel via 
simulations. For the given RFI scenario, several C programs have been written in 
order to simulate the convolutional coded bit error rate p for rate 1/2 convolutional 
code and its punctured higher rate cousins (e.g., rate 2/3, 3/4 and 4/5). As mentioned, 
if the RFI affects more than one channel code bit, the convolutional interleaver will 
be used against RFI. A brief introduction of this interleaving technique will be given 
later in this Section. 

The throughput simulation procedure for AWGN-plus-RFI channel is much like 
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that of the AWGN-only channel. The simulation is a combination of the both C 
programs and Matlab programs. Following is the throughput computation procedure 
for AWGN-plus-RFI channel. Again, we point out that Code 3 is assumed. 

Step 1 : Generate a time- varying Gaussian shaped E c (syN 0 channel as in Fig. 3.4. 
Step 2: For a sufficiently large number of E c /N 0 values (e.g., 21), determine p for 
the AWGN- plus-RFI scenario of interest via C program simulations, for 
code rates 1/2, 2/3, 3/4 and 4/5. Then produce p vs. symbol number (time) 
curves by interpolating these measured point. ( Repeat for other RFI 
scenarios.) 

Step 3 : Calculate P cw by (3.3), adjust the convolutional code rate and then the RS 
code rate so that P cw < 10‘ 5 . Store the composite rate r(s) for later use. 

Step 4: Calculate T a by (3.5), then compare with the theoretical throughput T K 
(e.g., 175,300 bits) to obtain throughput efficiency. 

Convolutional Interleaving Technique 

The idea of interleaving is very simple. Basically, the order of a sequence of 
symbols is rearranged in a deterministic manner before the sequence is sent to the 
burst-error channel, and then at the receiver side, the inverse permutation is applied to 
restore the sequence to its original ordering. 

The interleaving structure we shall discuss here was proposed by Ramsey [12]. An 
example of a convolutional interleaver and its corresponding de-interleaver with m = 
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4 delay lines and D = 1 symbol delay elements is shown in Fig. 4.1. As we can see, 
the input symbols are read onto the delay lines in a certain order, and the output of the 
delay lines is read in the same order. 

Because of the separation provided by the delay elements, any burst of b < m 
errors inserted by the channel results in single errors at the de-interleaver output. For 
example, in the figure, burst errors hit x 8 , x 5 , and x 2 , but those symbols are separated 
by "three at the output of de-interleaver. This implies that the good long burst-error- 
correction codes can be obtained from the good short burst-error-correction codes. 

With interleaving, extra delay is needed for the decoder. Defining the parameter 
L = mD (= m in our case in which Z>=1) 

the total end-to-end delay of convolutional interleaving is m(m- 1) symbols. Typically, 
the interleaver parameters m, D, and L are selected as follows 


m > the length of the burst error 
L > the decoding path length 

• • •, V X 3 , X . 6 , X 9 , X 4 , Xp X 2 , X 5 , Xjj, Xs, X 2 , x ,, x^, ... 



•••> Xp* X|> *2’ X 3’ X 4’ *5’ \< *T X 8’ V ^lO’ Nl’ ^2’ *0’ Xp ^2> Xj> X 4 , ^5> X^ *8, Xg Xj 0 , X, p X^, - • 


Fig. 4.1 Convolutional Interleaver and Corresponding De-interleaver 
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D = 1 


In this thesis, since the decoding path length is 32 for a rate 1/2 convolutional 
code and 64 for its higher rate (>l/2) punctured code, L has to be chosen larger than 
32. Since it is easier to use the same interleaver for all convolutional code rates in the 
implementation, we set the interleaver parameter m to m = 64 delay lines (although 
the rate 1/2 code really needs only 32 delay lines). 

4.2 Simulation Results 

In this section, the performance results of Code 3 for the AWGN-plus-RFI 
channel will be shown. Table 4.1 shows the RFI scenarios we examined in this thesis. 
In the table, MX is the mean interarrival time for the Poisson arrival model. We have 
Table 4.1 RFI Scenarios Examined for the Channel 


RFI (a = 0.1) 

RFI/code-bit ratio (3 

arrival model 

Case 1 

o.i 

Poisson (arrival rate X = 9.5e-3) 

Case 2 

0.1 

periodic (duty = 9.5 e-3) 

Case 3 

5 

Poisson (arrival rate X = 7.43e-5) 

Case 4 

5 

periodic (duty = 7.23e-5) 

Case 5 

10 

Poisson (arrival rate X = 4.24e-5) 

Case 6 

10 

periodic (duty = 3.9e-5) 

Case 7 

20 

Poisson (arrival rate X = 2.23e-5) 

Case 8 

20 

periodic (duty = 2.07e-5) 
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selected X in each P > 1 case so that p rate m (Poisson) « p rate m (Periodic) rather than 
X = duty cycle. 

Fig. 4.2 plots the Case 1 results analogous to Fig. 3.9 except we have included the 
p vs. symbol number curves based on simulations. The threshold of p (i.e., 2.39e-3) is 
the highest value allowable to achieve the P cw < 10" 5 requirement for the (255, 223) 
RS code, the most powerdul RS code. We note that p_4/5 is above this threshold line, 
so that rate 4/5 code is never selected as seen in the third plot. We also note that the 
p_l/2 curve becomes flat in the region of high E/N 0 , because, in the high E/N 0 
region, the RFI is more severe than the noise and dominates the BER p. Because the 
RFI/code bit duration ratio P is very small (i.e., P = 0.1), the effect of RFI is limited, 
and we obtain 86.68 percent of theoretical throughput (i.e., 175,300 bits) for this RFI 
scenario, compared to the 87 percent on AWGN-only channel. 

Fig. 4.3 plots the results for Case 2, which is the Case of periodic arrival RFI with 
P =0.1. All four rate of convolutional codes can be used in this Case, and we obtain 
86.13 percent of theoretical throughput this time. 

For the P = 5 Cases, Fig. 4.4 and Fig. 4.5 plot the results for Case 3 (Poisson 
arrival RFI) and Case 4 (Periodic arrival RFI), respectively. Only three rates (i.e., rate 
1/2, 2/3, 3/4) convolutional codes were simulated since only the two lowest rate codes 
of them can be used to achieve P cw < 10" 5 (see second and third plots in Fig. 4.4 and 
Fig. 4.5). Because RFI affects more than one code bit, interleaving is employed. For 
Case 3, we obtain 74.74 percent of theoretical throughput, and for Case 4, , we obtain 
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EcNo(dB) vs. symbol number for BPSK AWGN(+RFI) channel 



code rates and Ro bound vs. symbol number AWGIM+RFI channel 
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Fig. 4.2 Results for Case 1 
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77.48 percent of theoretical throughput. 

For the p = 10 Cases, Fig. 4.6 and Fig. 4.7 plot the results for Case 5 and Case 6, 
respectively. Again, only rate 1/2 and 2/3 convolutional codes can be used in Case 5 
(shown in Fig. 4.6), but rate 3/4 code can be also used in Case 6 (shown in Fig. 4.7). 
For Case 5, we obtain 73.64 percent of theoretical throughput, and for Case 6, we 
obtain 74.34 percent of theoretical throughput. 

" Finally, the Cases of P = 20 was examined, and Fig. 4.8 and Fig. 4.9 plot the 
results for Case 7 and Case 8, respectively. Again, only rate 1/2 and 2/3 convolutional 
codes can be used (shown in Fig. 4.8 and Fig. 4.9). For Case 7, we obtain 74.84 
percent of theoretical throughput, and for Case 8, we obtain 74.45 percent of 
theoretical throughput. 
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EcNo(dB) vs. symbol number for BPSK AWGN(+RFI) channel 
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EcNo(dB) vs. symbol number for BPSK AWGN(+RFI) channel 
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EcNo(dB) vs. symbol number for BPSK AWGN(+RFI) channel 
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Fig. 4.5 Results for Case 4 
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code rate and Ro bit error rate 


EcNo(dB) vs. symbol number for BPSK AWGN(+RFI) channel 
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code rates and Ro bound vs. symbol number AWGN+RFI channel 
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Fig. 4.6 Results for Case 5 
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Fig. 4.7 Results for Case 6 
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EcNo(dB) vs. symbol number for BPSK AWGN(+RFI) channel 
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EcNo(dB) vs. symbol number for BPSK AWGN(+RFI) channel 
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Fig. 4.9 Results for Case 8 
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Table 4.2 summarizes all the results discussed above and together with the coded 
daily throughputs. 

Table 4.2 Coded Throughput Efficiencies and Coded Daily Throughputs 


RFI 

Coded 

Daily Coded Throughputs 

Scenarios 

Efficiencies 

(Daily Theoretical Throughput = 1.5 Gbits) 

Case 1 

0.8668 

1.3 Gbits 

- Case 2 

0.8613 

1.29 Gbits 

Case 3 

0.7474 

1.12 Gbits 

Case 4 

0.7748 

1.16 Gbits 

Case 5 

0.7364 

1.1 Gbits 

Case 6 

0.7434 

1.12 Gbits 

Case 7 

0.7484 

1.12 Gbits 

Case 8 

0.7445 

1.12 Gbits 
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Chapter 5 


SUMMARY AND CONCLUSIONS 


The design of a reliable satellite communication link involving the data transfer 
from a small, low-orbit satellite to a ground station, through a geostationary satellite, 
was examined. In such a scenario, the received signal power to noise density ratio 
increases as the transmitting low-orbit satellite comes into view, and then decreases as 
it departs, resulting in a short-duration, time-varying communication link. The 
optimal values of the small satellite antenna beamwidth, signaling rate, modulation 
scheme and the theoretical link throughput (in bits per day) from the first sub-project 
were taken to check the practical coding schemes, both fixed rate and variable rate 
concatenated FEC coding schemes for the AWGN channel. The effect of RFI on the 
best coding scheme, which maximized the daily throughput while satisfying a 
prescribed probability of error requirement, among above schemes was examined. 

In Chapter 2, the discrete-time equivalent AWGN-plus-RFI model for simulation 
was described. In Chapter 3, the throughput efficiencies for the AWGN channel were 
numerically estimated. Three coding schemes, one fixed rate and two variable rate 
schemes were examined. It was found that the Code 3, which is the variable rate for 
both the RS codes and convolutional codes, could achieve 87 percent of the 
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theoretical throughput, equivalent to 1.3 Gbits/day based on the cutoff rate R 0 . Then 
in Chapter 4, the effect of RFI on Code 3 was examined. Eight different RFI 
scenarios, including both periodic arrival RFI and Poisson arrival RFI, RFI/code bit 
duration ratio value from 0.1 to 20 were considered. Interleaving technique was used 
to mitigate degradation due to RFI. It was found that the Code 3 could still achieve 74 
percent of the theoretical throughput, equivalent to 1.11 Gbits/day based on the cutoff 
rate R 0 , compared to 87 percent for the AWGN-only case. 

The final conclusion is easy to be made based on the above results. Code 3 is 
recommened for this time- varying satellite channel. Actually, the chip set for variable 
rate convolutional codes are available in the market. 
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APPENDIX A 


PUNCTURED CONVOLUTIONAL CODES 

In order to transfer more information data in the higher SNR region, higher rate 
(>l/2) convolutional codes were used in this thesis. But a straightforward application 
of Viterbi decoding to high-rate codes becomes very rapidly impractical as the coding 
rate increases. Fortunately, the high-rate codes can be obtained by the periodic 
elimination (i.e., puncturing) of specific code symbols from the output of a low-rate 
encoder [13]. The high-rate punctured code depends on both the low-rate code, called 
the original code, and on both the number and specific positions of the punctured 
symbols. The pattern of punctured symbols, described in a matrix form, is called the 
perforation pattern of the punctured code. 

In this thesis, the original code is the rate- 1/2, constraint length 7 convolutional 
code, and its punctured higher rates (i.e., rate 2/3, 3/4, 4/5) codes have been used. The 
perforation pattern matrixes for these codes are given [13] as follows, where 0 
indicates the punctured bits, 

D 1 

M 2 , , = 

2/3 |_1 0 

["110 
0 1 
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Fig. A.l and Fig. A.2 show an example of how to obtain a punctured rate-2/3 code 
with the perforation pattern matrix of M m from its original rate- 1/2, 4-state encoder. 
(The constraint length 7 code has too many states, 64, to demonstrate this principle). 

Fig. A.l gives the trellis diagram of a simple four states, rate- 1/2 convolutional 
code where every fourth symbol is punctured (indicated by X). By reading this trellis 
two branches at a time and redrawing it as in Fig. A.2, we see that it corresponds to a 



Fig. A.2 Trellis for Punctured Rate-2/3 Convolutional Code 
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4-state, punctured rate-2/3 code. 

The punctured code can be decoded by adding relatively simple circuitry to the 
codecs of the original low-rate 1/2 code. The added circuitry is used to insert dummy 
data (at the decoder) into the positions corresponding to the deleted code symbol 
based on their perforation pattern matrixes, where these dummy data are discarded by 
assigning them the same metric value (usually zero) regardless of the code symbol in 
the decoding process. The advantage of this decoding method is that Viterbi codecs 
for high-rate punctured codes is that the same rate 1/2 decoder may be used to decode 
the entire family of codes. 
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APPENDIX B 


MATLAB AND C PROGRAMS LIST 


In this section, we describe all the programs used in this thesis and then give a 
sample C simulation program for the AWGN-plus-RFI channel at the end. All the 
programs listed here can be found from a disk attached at the back of this thesis. 

B.l Programs List 

Matlab Calculation Programs 

t_cv_rs3.m This program computes the throughput efficiency for a fixed rate 1/2, 
constraint length 7 convolutional code concatenated with a (255, 223) 
RS code (i.e.. Code 1) on the AWGN-only time-varying channel. 
t_vr_cc3.m This program computes the throughput efficiency for the variable rate 
(1/2, 2/3, 3/4, 4/5), constraint length 7 convolutional codes 

concatenated with a (255, 223) RS code (i.e., Code 2) on the AWGN- 
only time-varying channel. 

vcr3.m This program computes the throughput efficiency for the variable rate 

(1/2, 2/3, 3/4, 4/5), constraint length 7 convolutional codes 

concatenated with a variable rate (255, k ) RS code (i.e.. Code 3) on 
the AWGN-only time- varying channel. 
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vcr3 RFI.m 


pb_vit.m 


binom.m 


mpskroli.m 


pcw.m 


This program computes the throughput efficiency for the Code 3 on 
the AWGN-plus-RFI time-varying channel. This program need to load 
the probability of error at the output of Viterbi decoder (i.e, p) data 
file generated by pb_vit.m. 

This program produces p vs. symbol number curves by interpolating 
these measured points, where the p data are generated by the C 
simulation programs for the variable rate convolutional codes on the 
AWGN-plus- RFI channel. 

This fuction program calculates the binomial coefficients, “n choose 
k”. It is used in the first four programs discussed above. 

This fuction program calculates the cutoff rate R 0 in units of 
information bits/channel use. The format of this function is : 
mpskrO(m, EsNodB). It is used in the first four programs above. 

This function returns the probability of codeword error for a (n, k) RS 
code over GF(q) when the channel symbol error probability is p. The 
format of this function is: Pcw(n, k, p). It is also used in the first four 
programs discussed above. 
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C Simulation Programs 


RFIl.cpp 


RFI2.cpp 


RFD.cpp 


RFI4.cpp 


This program simulates the probability of bit error for a rate 1/2, 
constraint length 7, soft decision viterbi decoder in a RFI 
environment. The Viterbi decoder uses the best metric with 32-bit 
memory, and the interleaving is not emploied. The output of this 
program will be written to a file, which will be then used in 
vcr3_rfi.m. 

This program simulates the probability of bit error for a rate 2/3, 
constraint length 7, soft decision punctured viterbi decoder in a RPI 
environment. The Viterbi decoder will take two stages from trellis - at a 
time, so that the decoding memory is 64-bit. The interleaving is not 
invoved. The output of this program will be written to a file too. 

This program simulates the probability of bit error for a rate 3/4, 
constraint length 7, soft decision punctured viterbi decoder in a RFI 
environment. The Viterbi decoder will take two stages from trellis at a 
time, so that the decoding memory is 64-bit. The interleaving is not 
invoved. The output of this program will be written to a file too. 

This program simulates the probability of bit error for a rate 4/5, 
constraint length 7, soft decision punctured viterbi decoder in a RFI 
environment. The Viterbi decoder will take two stages from trellis at a 
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time, so that the decoding memory is 64-bit. The interleaving is not 
invoved. The output of this program will be written to a file too. 
RFIlint.cpp This program simulates the probability of bit error for a rate 1/2, 
constraint length 7, soft decision viterbi decoder in a RFI 
environment. The Viterbi decoder uses the best metric with 32-bit 
memory, and the interleaving is emplyed. The output of this program 
will be written to a file, which will be then used in vcr3_rfi.m. 
RFI2_int.cpp This program simulates the probability of bit error for a rate 2/3, 
constraint length 7, soft decision punctured viterbi decoder in a RFI 
environment. The Viterbi decoder will take two stages from trellis at a 
time, so that the decoding memory is 64-bit. The interleaving is 
invoved. The output of this program will be written to a file too. 
RFI3_int.cpp This program simulates the probability of bit error for a rate 3/4, 
constraint length 7, soft decision punctured viterbi decoder in a RFI 
environment. The Viterbi decoder will take two stages from trellis at a 
time, so that the decoding memory is 64-bit. The interleaving is 
invoved. The output of this program will be written to a file too. 

B.2 Sample C Simulation Program 

For the AWGN-plus-RFI channel, Three different cases were considered during 
the programming. The first case is the Viterbi decoding without interleaving (RFI 
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affect less than one code bit) for both rate 1/2, constraint 7 convolutional code and its 
high-rate punctured codes, the second case is the Viterbi decoding with interleaving 
(RFI affect more than one channel code bit) for the original rate 1/2 convolutional 
code, and the last case, and also the most complicated one, is the Viterbi decoding 
with interleaving for the high-rate punctured codes. Because the first two cases can 
be considered as the sub-programs of the last case, we just give a sample, Viterbi 
decoding with interleaving for 3/4 punctured code, C simulation program here, and 
from this program, the reader may enable to either figure out or easily understand 
that of the first two cases. 


Sample C Simulation Program 

f-k-k-k'k'k'k’k'k'k'k'k-k'k'k'k'kle-k’k'k'k-k'kic'k’k-k'k'k-k-kic-k'k-k'k'k'k'k'kic'kic-k-k'k'k-k-k'k-k-k'k’k'k j 

/* Simulate a rate 3/4 soft decision punctured */ 

/* viterbi decoder reference to [13] . */ 

/* The definitions of parameters used in this */ 

/* program are either listed below or followed by */ 

/* the parameter at the first time being useed. */ 


/* */ 
/* state_bk -- shift register of convolutional encoder,*/ 

/* Ck coded bits, */ 

/* _Rk received bits, */ 

/* biterrs bit error counter, */ 

/* err_max maximum number of error to stop */ 
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/* the simulation, */ 

/* opt_n the index of the highest metrix, */ 

/* opt_met optimal (highest) matrix of trells */ 

/* decision equal to 1 when decoder */ 

/* making error, and 0 otherwise, */ 

/* Index_DL -- index of delay lines for both */ 

/* interleaver and de-interleaver, */ 

/* pattern perforation pattern matrix of rate */ 

/* 2/3 punctured code, */ 

/* data_up uper 32-bit of information bit stream, */ 

/* data_low -- lower 32-bit of information bit stream, */ 

/* path_up uper 32-bit surviving pathes of trellis,*/ 

/* path_low -- lower 3 2 -bit surviving pathes, */ 

/* dt_m extra delay buffer due to interleaver */ 

/* for information bits, */ 

/* nbits number of information bits have been */ 

/* sent through, .*/ 

/* Fi random phase of RFI, */ 

/* Tc code bit duration, */ 

/* beta RFI/Tc ratio, */ 

/* delta_W deferential frequency of RFI, */ 

/* Pcc probability of bit error of Viterbi */ 

/* decoder */ 

/* mean_cw mean of Poisson random variable, */ 

/* k_cw RFI interarrival time counter */ 

/* */ 
j * * * * * * * * *★************★******★**★*★*******'************* / 
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#include <stdio.h> 

#include <stdlib.h> 

#include <math.h> 

#include "c:\temp\data.h" //Contains Ec (s) N 0 _dB data and 7-bit 

parity lookup table 

#define POLYA 0133 //First generator polynomal of rate 

1/2, constraint 7 convolutional 
encoder 


# define 

POLYB 

0171 

//Second generator polynomal 

#def ine 

DELAY 

64 

//Decoding 

length of Viterbi 

#def ine 

TWOPI ' 

6 . 283185307 


#def ine 

M_INF 

-1 . 0e8 



#def ine 

D_LINE 

64 

//Number of 

delay lines for a 


convolutional interleaver 


long randws (void) ; //Random number generator 

double expo (double) ; //Interarrival time generator 

for Poisson arrival RFI 

void interleaver (unsigned, float *, unsigned); 

void awgn(float *, float, float *) ; //AWGN generator 

void deinterleaver (unsigned, float *, unsigned) 

void viterbi (float *, float *, unsigned long *, unsigned long 
*) ; // Viterbi decoder 

float T_int [D_LINE] [D_LINE] , T_de [D_LINE] [D_LINE] ; 

//Delay lines tables for interleaver and de- interleaver 
struct table { 

unsigned ps[128] ; 
unsigned long in[64]; 
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int outl [128] ; 
int out2 [128] ; 

} trellis; 

main ( ) 

{ 

unsigned state_bk, n, i, j, k, err_max, opt_n,- 
unsigned pattern [6] ={ 1 , 1 , 1 , 0 , 0 , 1} , decision, Index_DL; 
unsigned long data_up, data_low, path_up[64], nbits, bk, 
path_low [64] , dt_m [95] ; 

float EbNo_dB , EbNo, EcNo, Ec, No, stddev, m, Pcc[120], 

metric[64], biterrs, opt_met, Tau, Fi, Alfa, beta, Tc, 
low, up, delta_W, Ck[2], Rk[2]; 
double k_cw, mean_cw; 

FILE * f _W ; 

/* generate trellis lookup table */ 
for (n=0; n<64; n++) { //Current state 

trellis . in [n] = n >> 5; //Input bit 

for ( i=0 ; i<2 ; i++ ) { //Last bit of previous state 

trellis .ps [2*n+i] = ( (n&31) << 1) | i; //Previous state 

trellis . outl [2*n+i] = Partab [ ( (n<<l) | I) & POLYA] ? 

1 : -1; //First output symnol 
trellis . out2 [2*n+i] = Partab [( (n<<l) | i) & POLYB] ? 

1 : -1; //Second output symbol 

} 

} 

Tc~~= l./4.5e6; //Because the data rate = 4.5e6 cbps 
Alfa = .1; // Set Alfa = 0.1 
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delta_W = Alfa*TWOPI/Tc ; // Delta omiga of CW RFI 
mean_cw = 23567.; k_cw = expo ( 1 . /mean_cw) ; 

Tau = 10*Tc ; beta = Tau/Tc; err_max = 80; 

printf("\n\n Program -is running, please wait 

\n \nBit Error Rate For Parameter Group A :\n\n"); 

for ( j = 0 ; j < 1 1 ; j++) { //ll sample E c /N 0 values in region 

[-15,000, 0] 

/* initialize some parameters */ 

state_bk = biterrs = nbits = 0; 

Index_DL = data_up = data_low = 0 ; 
for (i=0; i<95; i++) dt_m[i] = 0; 

//Need 2*63*3/4 32-bit input data buffer 
for (i=0 ; i<D_LINE ; i++) { 

for (m=0; m<D_LINE; m++) { //Preset 1 
T_int[i] [m] = T_de[i] [m] = 1; 

} / /For m loop 
} //For i loop 
m = 0 ; 

/* calculate noise level */ 

EcNo = pow (10., EcNo_dB [ j ] / 1 0 . ) ; 

Ec = 1 . ; 

No = Ec/EcNo ; / / AWGN level 

stddev = sqrt (No/2 . ) ; //Standard deviation of AWGN 
/* initialize metrics*/ 

for (n=0; n<64; n++) metric [n] = M_INF ; 

~ metric [0] = M_INF/100; 
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if ( (k_cw<2 ) && (k_cw>=l) && (k ! =1) ) { 

//The 1st RFI hit Rk[l] 

Fi = TWOPI* (randws () +1) /( (float) 0x100001) ; 
low = k_cw*Tc; //Lower limit of integrator 
if ( (k_cw+beta) < (m+1) ) { 

up = (k_cw+beta) *Tc ; //Upper limit 
Rk[l] += 10/ (Tc*delta_W) ) * (sin (delta_W*up+Fi) 
sin (delta_W*low+Fi) ) ; 

k_cw = expo ( 1 . /mean_cw) + up - m*Tc; 
m = 0; // Reset m=0, and exit the loop 
} // For if ((k_cw<2) && (k_cw>=l) && (k ! = 1 ) ) 

else { 

up = (m+1 ) *Tc ; 

Rk [1] += 10/ (Tc*delta_W) ) * (sin (delta_W*up+Fi) 
sin (delta_W*low+Fi) ) ; 

m++ ; 

} 

k_OW-- ; 

} 

else if (k_cw < 1) { 

//The 1st RFI hit Rk[0] or continue from above 
if (m == 0) 

Fi = TWOPI* (randws () +1) /( (float) 0x100001) ; 
low = (k_cw > m) ? k_cw*Tc : m*Tc; 
if ( ( (k_cw+beta) < (m+1) ) && ( k ! =2 ) ) { 

up = (k_cw+beta) *Tc ; 

Rk[0] += 10/ (Tc*delta_W) ) * (sin (delta_W*up+Fi) 
sin (delta_W*low+Fi) ) ; 

k_cw = expo (1 . /mean_cw) + up - m*Tc; 
m = 0 ; 


64 



} 

else if (k!=2) { 

up = (m+1 ) *Tc ; 

Rk [0] += (10/ (Tc*delta_W) ) * (sin (delta_W*up+Fi) 

- sin (delta_W*low+Fi) ) ; 

m++ ; 

low = (k_cw > m) ? k_cw*Tc : m*Tc; 
if ( ( (k_cw+beta) < (m+1) ) && (k!=l)) { 

up = (k_cw+beta) *Tc ; 

Rk[l] += 10/ (Tc*delta_W) ) * (sin (delta_W*up+Fi) 

- sin (delta_W*low+Fi) ) ; 

k_cw = expo ( 1 . /mean_cw) + up - m*Tc; 
m = 0 ; 

} 

else if (k!=l) { 

up = (m+l)*Tc; 

Rk [ 1 ] += 10/ (Tc*delta_W) ) * (sin (delta_W*up+Fi) 

- sin (delta_W*low+Fi ) ) ; 

m++ ; 

} 

. } 

} 

else k_cw -= 2; 

/* deinterleaving the received bits */ 
deinterleaver (Index_DL, Rk, k) 

Index_DL = k ? ( Index_DL+l) %D_LINE : 

( Index_DL+2 ) %D_LINE ; 

/* viterbi decoder */ 

if (nbits > D_LINE* (D_LINE-1) *3/4) { 
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Rk [0] *= pattern [2*k] ; // Insert dummy data 

Rk[l] *= pattern [2*k+l] ; 

viterbi (Rk, metric, path_up, path_low) ; 


/* check errors */ 

if (nbits >= int (DELAY+D_LINE* (D_LINE-1) ) ) { 

opt_met = M_INF; 
for (n=0 ; n<64; n++) { 

if (metric [n] > opt_met) { 
opt_met = metric [n] 
opt_n = n; 

} 

} 

decision = (path_low [opt_n] A data_low) & 1; 
biterrs += decision; 

} 

} 

} 

Pcc[j] = biterrs/(nbits-int(DELAY+D_LINE*(D_LINE-l)*3/4)); 
printf ( "%10 . 3e, ", Pcc[j]); 

/* open a file for write */ 

if((f_w = fopenCp3_c4.dat ","a")) == NULL) { 

printf ("\n The p3_c4.dat file cannot be opened\n"); 
return 0 ; 

} 

f print f(f_w, "%10 . 3e, ", Pcc[j]) ; 

/* close the file */ 
f close (f w) ; 
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} 

print f ( " \n\n\nThe End ! " ) ; 
return 0 ; 

} 

/* Random # generator described by Widrow and Stearns */ 
long randwsO 
{ 

static long seed = 78971; //seed to arbitrary value 
return! seed = (2045*seed+l) & Oxfffff ); 

} 


/* Exponential random variable generator described by ee572 
text book */ 

double expo (double a) //generate the interarrival time of 
Poisson arrival RFI 
{ 

return ( -log ( (randws ( ) +1) / ( (float) 0x1000 01) ) /a) ; 

} 

/* AWGN generator */ 

void awgn(float Ck[], float stddev, float Rk[]) 

{ 

float ul, u2, a, b; 

ul = (float) (randws () +1) /( (float) 0x100001) ; 
u2 = (float) (randws () +1) /( (float) 0x100001) 
a=TWOPI*ul ; 

b^sqrt ( -2.0*log(u2) ); 

Rk [0] = Ck[0] + stddev*b*cos (a) ; 
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Rk[l] = Ck[l] + stddev*b*sin (a) ; 

} 

/* Viterbi decoder */ 

void viterbi ( float Rk[], float metric[], unsigned long 
path_up[], unsigned long path_low[]) 

{ 

unsigned n, i; 

unsigned long path_up_new [64] , path_low_new [64] ; 
float m[64] [2] ; 

for ( n= 0 ; n< 6 4 ; n++) { // Current state 

for ( i = 0 ; i<2 ; i++) // Input bit for previous state 

m [n] [i] = metric [trellis .ps [2*n+i] ] + 

Rk [0] *trellis . outl [2*n+i] + Rk [1] *trellis . out2 [2*n+i] ; 
if (m [n] [0] > m [n] [1] ) { 

path_up_new [n] = (path_up [trellis .ps [2 *n] ] >> 1) | 

(trellis . in [n] << 31); 

path_low_new [n] = (path_low [trellis .ps [2*n] ] >> 1) 

| ( (path_up [trellis .ps [2*n] ] & 1) << 31); 

} 

else { 

path_up_new [n] = (path_up [trellis .ps [2*n+l] ] >> 1) 

| (trellis . in [n] << 31); 

path_low_new [n] = (path_low [trellis .ps [2*n+l] ] >> 

1) | ( (path_up [trellis .ps [2*n+l] ] & 1) << 31); 

} 

} 

for (n=0; n<64; n++) { 

metric [n] = (m[n] [0] > m[n] [1]) ? m[n] [0] : m [n] [1] ; 

// Metric update 
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path_up [n] = path_up_new[n] ; // Path_up update 

path_low[n] = path_low_new [n] ; // Path_low update 

} 

} 

} 

/* interleaver */ 

void interleaver {unsigned Index, float Ck[], unsigned k) 

{ 

unsigned n; 
float temp; 
if (k ! = 2) { 

temp = Ck[0]; // update CK[0] 

if (Index >0) { // update delay line of the table 

for (n=Index; n>0; n--) 

T_int [Index] [n] = T_int [Index] [n-1] ; 

} 

T_int [Index] [0] = temp; 

Ck[0] = T_int [Index] [Index] ; 

} 

if (k != 1) { 

if (k==0) Index = (Index+1) %D_LINE; 

// (index+1) maybe back to 0 
temp = CkTl]; // update CK[1] 
if (Index >0) { 

for (n=Index; n>0; n--) 

T_int [Index] [n] = T_int [Index] [n-1] 

} 

~ T_int [Index] [0] = temp; 

Ck[l] = T int [Index] [Index] ; 
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} 


} 

/* deinterleaver */ 

void deinterleaver (unsigned Index, float Rk[], unsigned k) 

{ 

unsigned n; 
float temp; 
if (k != 2) { 

temp = Rk [ 0 ] ; // update RK[0] 

if ( (D_LINE- Index- 1) > 0 ) { 

for (n= (D_LINE- Index-1) ; n>0; n--) 

T_de [Index] [n] = T_de [Index] [n-1] ; 

} 

T_de [Index] [0] = temp; 

Rk [0] = T_de [Index] [D_LINE- Index- 1] 

} 

if (k ! = l) { 

if (k==0) Index = (Index+1) %D_LINE; 

// (index+1) maybe back to 0 
temp = Rk [ 1 ] ; // update RK[1] 
if ( (D_LINE- Index- 1) > 0 ) { 

for (n= (D_LINE- Index-1) ; n>0; n--) 

T_de [Index] [n] = T_de [Index] [n-1]; 

} 

T_de [Index] [0] = temp; 

Rk [1] = T_de [Index] [D_LINE- Index- 1] ; 

} 

} ~ 

^•k-k-k'k-k’k-k-k'k'k'k'k'k'k Elld Of P2T0^3r3.m ****************** j 
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